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The rates governing momentum broadening in a quark-gluon plasma with a mo- 
mentum anisotropy are calculated to leading-log order for a heavy quark using kinetic 
theory. It is shown how the problematic singularity for these rates at leading-oder 
is lifted by next-to-leading order gluon self-energy corrections to give a finite con- 
■ tribution to the leading- log result. The resulting rates are shown to lead to larger 

. momentum broadening along the beam axis than in the transverse plane, which 

PsJ \ is consistent with recent STAR results. This might indicate that the quark-gluon- 

£^ ' plasma at RHIC is not in equilibrium. 

I. INTRODUCTION 

One of the hottest debated questions in the context of ultrarelativistic heavy ion collisions 
^ ■ is whether a "thermal" quark-gluon plasma has been created at the highest energy runs of the 
Relativistic Heavy- Ion Collider (RHIC). On the one hand, the success of ideal hydrodynamic 
fits to experimental data have been interpreted [l[ to imply rapid thermalization of the bulk 
q \ matter at RHIC. On the other hand, perturbative estimates of the thermalization time 
result in much larger values, seemingly excluding the possibility of rapid thermalization (see 
however |3j) . Plasma instabilities may help speeding up the equilibration process b y rap idly 
isotropizing the system [3, H, @, 0, [1] , probably by some cascade to the UV [§, [l(J HI, H, 
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although they might still be too slow to overcome the initially strong longitudinal expansion 
rate |l4|. Thus, it is currently unclear how (and if) the quark gluon plasma created at RHIC 
stays locally isotropic (let alone thermal) during most of its time evolution. 

Nevertheless, many theoretical calculations in the context of heavy-ion collisions assume 
a homogeneous and locally isotropic system. However, the physics of anisotropic plasmas 
(even if they are close to isotropy) can be quite different from that of isotropic plasmas, c.f. 
the presence of plasma instabilities in the former. Due to these potentially large differences, 
it might be interesting to reanalyze calculations of experimentally accessible observables in 
the context of anisotropic plasmas, wherever possible. 

In what follows, I will thus study momentum broadening in a homogeneous but locally 
anisotropic system. Since I will be interested in qualitative effects only, I limit myself to con- 
sidering momentum broadening for a heavy quark induced by collisions to leading logarithm 
(LL) accuracy. In section [Til I will discuss the calculational setup for anisotropic plasmas. 
Section UTTl will contain the calculation of the rates governing momentum broadening, which 
can be done analytically in the small anisotropy limit. In section ITVl finally I will try to apply 
these results to preliminary STAR data on jet shapes and discuss possible consequences. 



II. SETUP 

The setup of my calculation will be based on the work by Moore and Teaney [3]. They 
considered a heavy quark with momentum p moving in a static (thermal) medium with a 
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temperature T and I will quickly repeat some equations from their work which I want to 
use in the following. 

Due to collisions, the quark will loose momentum and the variance of its associated 
momentum distribution will broaden according to 15 
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(1) 



where I choose my coordinate system such that ((Ap\\) 2 ) = (p|| — (p\\)) 2 is the variance 
of the momentum distribution in the direction parallel to the direction of the quark and 
((Ap_i_) 2 ), ((Ap z ) 2 ) are the variances transverse to the direction of the quark (see also FigJT]). 

The functions tjd, k\\, k±, k z which encode average momentum loss as well as transverse 
and longitudinal fluctuations are calculated using kinetic theory. Schematically, they are 
given as [15j 
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where v = p/p°, f(k) is the gluon (quark) distribution function, M. is the scattering matrix 
and f k = J (similarly for q). If the transferred energy q° ~ v ■ q is small, one can 
approximate 



f(k) l±f(k-q°) - f(k - q°) [1 ± f(k)] ~ -%f{k)[l±f{k)] 



f(k) l±f(k-q c 



T ' 

f(k) [1± /(*)]. 



(3) 



Thus, if the quark is non-relativistic or if q° T, the coefficient km can be related to the 
energy loss rate — ^ by K|| = — f?^- This implies that also the other coefficient functions 
can be calculated by the appropriate modifications in the integrand of the energy loss rate. 



A. Collisional Energy Loss 



Studies of collisional energy loss of a heavy quark in isotropic systems have a long history 



[H, H E3, m, 0, 0, H 0H0. Assuming a perturbative expansion in powers (and 
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logarithms) of the strong coupling a s — |-, the leading-order logarithmic contribution to 
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the collisional energy loss can be obtained by calculating the matrix elements in Eqns. (J2]) 
with HTL propagators and restricting to soft momentum transfer q° < T where Eqns. (J3]) 
can be used [15 . 

There is, however, an alternative way by Thoma and Gyulassy [l8[ to calculate the 
collisional energy loss which is equivalent to calculating the HTL matrix elements and doing 
the integrals in Eqns. (CD). It is based on the expression for the energy loss of a quark by 
interaction with its induced chromoelectric field, 



dp 
~~dt 



Re|d 3 *J^(*,x).E? nd (*,x). (4) 

Here J" xt = g a v5 3 (x — vt) is the current of the quark with color charge q a and velocity v 
and .Kind is the induced electric field. In Fourier space Q = (q, to), E^d can be written as 

El nd (uj, q) = iu (GV - G#) JUQ), (5) 

where G^ and Gq are the full and free retarded gluon propagator, respectively. 
The collisional energy loss for a theory with N c colors thus becomes 

where the momentum exchange is restricted to |q| < T and only the leading logarithm should 
be kept. Because this restriction corresponds to Eqns. ([3]), the fluctuation coefficients may 
be glossed directly from Eqns. ([1]), finding 

# 2 (A 2 -1) T f d 4 Q 2Tqf { , , 
S 2 (A 2 -1) T f d 4 Q 2Tg 2 



-Im / 44^^ {G ij - vi(2n)5(u - v • q). (7) 



A convenient feature of these Equations is that all the information about the medium 
resides in the gluon propagator G. Thus, in this form Eqns. ([7]) are valid both for isotropic 
as well as for anisotropic systems. 



B. Anisotropic Plasmas 

In the rest frame of a thermal system, particles move in all directions with equal prob- 
ability. However, in the rest frame of a non-equilibrium system, particles might e.g. move 
predominantly in a plane, having an anisotropic probability distribution. Thus, dropping 
the assumption of isotropy, there will be at least one preferred direction in the system, which 
in the following I will take to be the z-direction. In the context of heavy-ion collisions, one 
can identify this direction with the beam-axis along which the system expands initially. As 
a model for an anisotropic system, I will assume that the anisotropic quark and gluon dis- 
tribution functions /(k) are related to the usual (isotropic) Fermi and Bose distributions 
/;>„( k ) as [3 

/(k) = N(Of iso (^k 2 + £(k ■ e.)2) , (8) 
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where e z denotes the unit- vector in the z-direction, £ is a parameter controlling the strength 
of the anisotropy and iV(£) is a normalization factor which in the following I will set to one. 
Note that £ = corresponds to the case of an isotropic system. 

An anisotropic system by definition does not have one temperature. However, it may have 
a dimensionful scale that is related to the mean particle momentum and which takes over 
the role of a "hard" scale in much the same way the temperature does in a thermal system. 
For a quark gluon plasma (even out of equilibrium), the natural choice seems to be the 
saturation scale Q s . Since kinetic theory calculations remain valid even out of equilibrium 
as long as there is a separation of scales between hard and soft momentum modes, one could 
replace T by Q s times some coefficient when translating formulae from an isotropic to an 
anisotropic system. For simplicity, however, I will refrain from doing that and keep T as a 
placeholder of the correct hard scale even for anisotropic systems. 

As outlined before, the main difference between calculating the fluctuation coefficients 
Eqns. ((7j) in isotropic and anisotropic systems, respectively, is the form of the gluon prop- 
agator G. While for isotropic systems the propagator in the HTL approximation is given 
by 1 
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one finds for systems with an anisotropy of the form Eq. (jSJ) within the same approximation 
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G ij (Q) 



+ 



(q 2 -uj 2 + a(u, q) + 7(0?, q)) q l q j /q 2 



q 2 - uj 2 + a(uj, q) ' (q 2 - uj 2 + a(u, q) + 7(0?, q)) (f3(uj, q) - uj 2 ) - q 2 n 2 <5 2 (u;, q) 

(P(u>, q) — oj 2 ) n*n J /n 2 — 5(u, q) (q l n j + q j n l ) 
(q 2 - uj 2 + a(u, q) + j(uj, q)) (f3(uj, q) - uj 2 ) - q 2 n 2 5 2 (co>, q) ' 



where here n = e z — qq^/q 2 and the structure functions a, (3, 7, 8 are generalizations of 
and to anisotropic systems. In the case of small anisotropy £, they are given as [261 ] 



a 



n T (c£;)+£ 
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'•*J /_ «9 1N 9 I a9 2 
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(3g 2 - l))/^ 2 , 



+n L (^) ((2g 2 -l)-^ 2 (3g 2 -r 



1 Here and in the following (unless stated otherwise) I use the specific choice of gauge A = 0; physical 
observables are not affected by this choice as they are manifestly gauge- invariant. 
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7 = |(3n r (^)-m^ 2 -l)(l-g z 2 ), 

S = ^LJ 2 m 2 D + m T {u>){l - Au 2 ))q z , (11) 

where u = uj/\q\ and q 7 , = q z /\q\. For larger values of £ one has to resort to numerical 
evaluations of a, j3, 7, 5 [261 ] . 

In all cases is the isotropic Debye mass, 

Q 2 f°° 

m 2 D =^—\ dppfUlPl), (12) 
7T Jo 

which for a thermal system with N c colors and Nf light quark flavors becomes 

»4 = *5*(i + ¥). 



III. MOMENTUM BROADENING 

As can be quickly verified, the general form of the propagator Eq. (flOl) reduces to the 
simple form Eq. in the limit of vanishing anisotropy £ —>■ 0. Thus, regardless whether 
the system is isotropic or not, the fluctuation coefficients are found by inserting Eq. (JTUJ) 
into Eqns. ([7]). 

The frequency integration is trivial, whereas the integration over total momentum ex- 
change |q| contains a subtlety for non-vanishing anisotropy £: since the static limit (a; — > 0) 
of e.g. a is real and negative, the propagator develops a singularity for space-like momenta 
Indeed, this singularity signals the presence of instabilities in the system (25I, 26[. Since 



|q 

these instabilities correspond to soft gauge modes that grow nearly exponentially until they 
reach non-perturbative occupation numbers, one might dismiss any perturbative approach 
such as mine as futile. However, my calculation should be applicable up to the time where 
the soft mode occupation number has not grown non-perturbatively large yet. Even with 
fast growing instabilities, this can be a long physical timescale if either the initial fluctua- 
tions are tiny or there is a counter-acting effect such as the expansion of the system. In the 
latter case, one can even argue that unstable modes do not grow non-perturbatively large 



during plasma lifetimes at RHIC |14|]. Therefore, even though for quantitative results one 
probably has to resort to numerical simulations, ignoring the effects of the non-perturbative 
soft gauge modes may be not such a bad approximation at least for some period of time. 

As a consequence, here the propagator singularity is mainly a practical obstacle to cal- 
culating physical observables since this singularity is in general non-integrable, as was first 



recognized by Arnold, Moore and Yaffe [25 



However, the fact that a has an imaginary part that vanishes only linearly in the static 
limit cures this problem for certain observables (namely those with numerators that also 



turn out to vanish in the static limit), which has been dubbed "Dynamical Shielding" [27 

In the formulation I have chosen, Eq. ([7]), there is only one power of the propagator and 
consequently the singularity is always integrable. However, the problem resurfaces in the 
next integration step, where for some of the fluctuation coefficients there is an uncanceled 
u in the denominator of the integral. Unlike the collisional energy loss, which is finite 
to full leading-order (LO), dynamical shielding breaks for at least some of the fluctuation 
coefficients. I will show in the following how to calculate these to LL accuracy. 
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FIG. 1: Sketch of the coordinate system: the quark moves along the || direction, perpendicular to 
the beam axis which coincides with the anisotropy direction. 

The calculation is divided up into two parts: first I calculate the "regular" contribution 
within the LL approximation, which involves only quantities for which dynamical shielding 
works. Then I will deal with the contribution which would give the LO contribution (the 
constant under the log) in an isotropic system and show that in anisotropic systems, it gives 
a contribution to the LL. 



A. Regular contributions 



Limiting oneself to LL approximation, the |q| integration in Eq. (J7j) can be done as in 
Ref. 27|, and one finds e.g. 
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(13) 



where q± = q±/\q\, 5 = S/\q\, Cj = v ■ q/|q|. 

I will consider the situation where the quark velocity is perpendicular to the beam axis, 
v • e z = 0. Thus I have Cj = vq\\, n ■ v = ujq z and n 2 = 1 — q 2 . From the explicit form of 
the small £ structure functions a, /3, 7, S in Eq. (fTTj) it can be seen that their imaginary part 
always involves at least one power of Cj. All the Cj's in the denominator of the integrand in 
Eq. (fT3"j) are thus canceled, showing indeed that there is no singularity to LO accuracy. 

For small anisotropies, the remaining integrations can be done analytically using 
Eqns. (MTT|) . I find 
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(14) 



As expected, these expression reduce to the known results [15[ in the isotropic limit (£ — > 0). 



B. Anomalous Contribution 

Let us now concentrate on the naive "constant under the log" contribution of Eqns. (1141) . 
To this end, it is illustrative pick out the first term of the propagator Eq. ( flOl) . and pluck it 
into Eqns. ([7]), finding 

(2tt) 3 ^V Hz I q 2 -cj 2 + a(w,q) V 7 

where the index % denotes i =_L, 1 1, z and I have used the same notations as in the previous 
subsection. Denoting this part of the propagator as A^ 1 = q 2 — uj 2 + a(cu, q) allows me to 
rewrite 

«i ~ / -0^Q"m 2 D f(q z ) (v 2 -Cj 2 - q 2 z Cj 2 ) qA A A* A , (16) 

where I introduced mjjCjflqz) = Ima and the * means complex conjugation. This form now 
directly involves the squared matrix from Eqns. (T5]). As discussed before, the static limit of 
a is real and negative to leading order in g, which produces a non-integrable singularity in 
Eq. f[T6"j) (unless i — || for which q\\v = Cj cuts off the singularity). 

However, from the structure of Eq. f[T6"j) . this singularity is reminiscent of so-called pinch- 
ing singularities, which are usually due to incomplete resummations of the propagator. In- 
deed, I find it is plausible that a has a non-vanishing imaginary part in the static limit at 
order 2 0(g 3 ), which allows one to integrate 
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dqq 3 \A A \ 2 ' 



2(1 -Cj 



2\2 



Re a . /Re a + g 2 (l - Co 2 ) 

■Arctan 



Ima y Ima J 
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Cj 2 ) 2 



(17) 



The second term in this equation gives a contribution to the LL for q = T ^> nip which 
has been accounted for already in the previous subsection. The other term, however, now 
involves Ima rather than the Cj of Eq. ( fl5l) in the denominator, which because of the 0(g 3 ) 
contribution behaves as 

lim Ima~ m 2 D (Cj + eg) , (18) 



2 See the appendix for details. 
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where c is just a number. Thus, the singularity at Cj — > 0+ is cut-off by this higher-order 
resummation in the self-energy. However, this also entails that the naive LO correction gives 
a contribution to the LL instead, as I will show in the following. 

First note that for any nonzero Cj, the first term on the r.h.s. of Eq. (1171) is finite also for 
c = and thus only contributes to LO. In order to extract a potential LL contribution, one 
may thus take Cj — > + everywhere except for the Cj in the denominator, which one has to 
replace by Cj — ► Cj + c g. For simplicity, I do this already in Eq. (1151) . using 3 



lim ImA^ 



-7t5 (q 2 + Re 



a 



(19) 



As can be verified, for small anisotropies £ the leading contribution to K{ is entirely given 
by this part of the propagator 4 . The integral over Cj then gives log (g) together with some 
finite LO contribution that I ignore. 

Thus, upon identifying log (g) = — log I find for the anomalous contributions to ft. 
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such that together with the regular contributions Eq. (JHJ), the full LL fluctuation coefficients 
in the small anisotropy limit are given by 



K± 



9 2 {Nl-l)Tml log T 



2N r 



An mo 



2Av 5 



(-3u + 8? 



1 (l-^ 2 )\ 
— + ; Arctanh {v) 

2v A 2v 6 



13w 5 - 2v 6 + 3(1 - v 2 ) 3 Arctanh(w 



2N r 



I) Tm 2 D , T 
— log 



47T 

24w 5 



1 

2^2 



(1 



,2\2 



2w 3 



-Arctanh (f) 



-9v + 24w 3 - 7w 5 - 6w 6 + 9(1 



,2\3 



Arctanh(f 



(21) 



C. Larger Anisotropies 

The relevant integrals of the previous subsections may also be evaluated numerically for 



arbitrary £ using techniques from 26|, |27| . The results are shown in Table [H For sufficiently 
small £, the results turn out to coincide with the analytic results Eqns. (|2"T|) . 

From Table [J it can be seen that the anomalous contribution makes up only a few 
percent at small £, while it becomes more important for larger anisotropies. Indeed, at 
larger £, besides Eq. ( I19p then also the other parts of the propagator contribute, further 



3 This (5-function is not to be confounded with the structure function 5\ 

4 This is because the structure function S 2 in the denominator may be dropped to leading order in 0(£). 
However, this result is specific to choosing the quark velocity perpendicular to e z , so that n • v — > in 
the static limit. 
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TABLE I: The coefficients multiplying y v 2 £ — - 4?r log in n±, k z for larger values of the system 
anisotropy £, evaluated numerically. For convenience, also the regular contributions are given 
explicitly. 



enhancing the anomalous contributions. For the total k Z} where the anomalous contribution 
is stronger, this manifests itself by a decrease as a function of v for larger £, whereas the 
regular contribution would have had the inverse trend. 

Another trend that is apparent from Tabled] is that the ratio k z / k± for a given anisotropy 
£ decreases for increasing velocity v. Put differently, the higher the quark's momentum, the 
more circular its associated jet shape. 



IV. DISCUSSION 

Both from the analytic results Eqns. (12 ip as well from the numeric evaluation given in 
Table [U one finds that for anisotropic systems k z /k± is always larger than one. In other 
words, a charm quark jet in an anisotropic quark-gluon plasma will generically experience 
more broadening along the longitudinal direction than in the reaction plane. Assuming an 
initial jet profile with ((Ap±)l) = ((Ap z )l), the ratio k z /k± can roughly be associated with 
the ratio of jet correlation widths in azimuth (A0) and rapidity (A77) as 

- gj. (22) 

Therefore, the calculation in the previous sections generically predicts this ratio to be larger 
than one in anisotropic systems. 

Due to the approximations in my calculation (heavy quark, leading-log, only collisional 
broadening), many effects that will change both k z and k± for real jets have been ignored. 
However, inclusion of other effects should not change my result on the ratio of k z and k± 
qualitatively, since k z = k± in isotropic systems. Thus, k z /k± > 1 in anisotropic systems 
should hold true also for real jets. Interestingly, recent measurements of the dihadron 



correlation functions by STAR [28] for Au+Au collisions at y's = 200 GeV indeed seem to 



imply much broader correlations in A77 than in A0. In FigJ2], STAR data for these correlation 




FIG. 2: Dihadron correlation functions in azimuth A(j) and space-time rapidity At/, for d+Au 
and central Au+Au at \fs = 200 GeV (figure courtesy J. Putschke (STAR), Proceedings of Hard 
Probes 2006). For comparison, on the right plot I show an ellipse with eccentricity e ~ which 
corresponds to a ratio k z /k± ~ 3. 



is compared to a general ellipsoidal shape with eccentricity e ~ which according to 
Eq. (122j) I would associate with a ratio k z /^± ~ 3. 

In other words, a charm quark jet with a "mean" momentum p ~ 4.6GeV (or v ~ 0.95) 
and a "mean" system anisotropy £ ~ 10 should experience broadening roughly consistent 
with the ellipsoidal shape in FigfSJ where "mean" here is referring to a mean over the whole 
system evolution. To get more quantitative, one would have to calculate the momentum 
broadening dynamically accordin g to Eqns. (pQ), which can e.g. be done within a viscous 
hydrodynamics simulation 29, 30, 3l[. However, note that as a general feature of anisotropic 



plasmas, one seemingly obtains jet shapes that can be much broader in rapidity than in 
azimuth. Other effects distorting the shape of jets in an isotropic plasma, namely flow effects 
(32| , probably cannot account for such dramatic asymmetries unless invoking extreme flows. 

Interestingly, momentum-binned data of the dihadron correlations from STAR seem to 
show a trend towards more circular jet shapes for higher trigger momentum, which I also 
found in section UTTlC. Precise data could thus provide further tests for calculations, maybe 
helping to constrain the system anisotropy. 

In view of this, a finite system anisotropy could be a natural explanation for these broad 
rapidity correlations in central Au+Au collisions at \/s = 200 GeV. Thus, turning the 
argument around, it may be that the sizable ratio (At/) / ( A0) is an indication that the 
plasma created at RHIC is after all not in equilibrium during a sizable fraction of its lifetime, 
calling into question the validity of the strongly advocated "perfect fluid" picture. 

Clearly, in order to shed more light onto this issue, many caveats on the application of my 
result to experimental data have to be addressed. Among other things, a calculation of the 
full LO correction as well as results for light quarks, gluons and the inclusion of radiation 
effects in anisotropic plasmas would be on the wish-list. Nevertheless, I hope that I was able 
to highlight the feasibility and potential value of re- calculating experimentally interesting 
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FIG. 3: Diagrams contributing to the gluon self-energy 



observables for anisotropic plasmas. 

To summarize, I have calculated the rates governing momentum broadening for a heavy 
quark to leading-log accuracy in an anisotropic quark-gluon plasma. The potential singu- 
larity, naively plaguing the constant under the log of these rates is shown to be cured after 
resummations of the NLO gluon self-energy, eventually contributing to the leading log. I 
expect this procedure to render finite also other observables that are affected by these sort of 
singularities, making their calculation feasible also in anisotropic plasmas. The rates found 
in this calculation suggest jet shapes that are elongated along the rapidity direction, which 
could explain recent STAR data and thus point towards the possibility of non-equilibrium 
phenomena at RHIC. 
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APPENDIX A: NLO GLUON SELF-ENERGY 



I will be interested in the static limit of the imaginary part of the gluon self-energy to 
0(g 3 ). Following Braaten and Pisarski [33J, this can be calculated by dressing the propaga- 
tors and vertices of the contributing diagrams. 



A further simplification occurs due to the static limit: following Rebhan [34J, only static 
loop momenta have to be resummed and no fermionic loops contribute. Since the ghost 
self-energy vanishes at leading order, the only contribution to the gluon self-energy at 0(g 3 ) 
in the static limit are thus the two diagrams shown in Fig. [3J 

Of these, the tadpole diagram (shown in Fig 121(a)) with bare vertices contributes [3 l 



SUT = -g 2 N c 



d 3 k , 
-, — tt7 lim 

(2tt) 3 



<T(G>,k) - Gg )o (u;,k)) - (G">,k) - G%>,k)) 



(Al) 

where G and G denote the dressed and free propagator, respectively. The dressed propa- 
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gator in covariant gauges can be found by inverting the relation 

9^K 2 + (-- l)K„K v + IV = G~l, (A2) 

where p is the gauge parameter and LP" is the gluon self-energy to 0(g 2 ) which for isotropic 
as well as anisotropic systems can be expressed [26| by the functions a, /3, 7, 5 introduced in 
the main text. 

For isotropic systems, inversion of Eq. (lA2j) is trivial and one indeed recovers the known 



result [35[ for the tadpole contribution to 511. 

For anisotropic systems, the inversion of Eq. (1A2|) is a little more involved but it is 
readily performed by e.g. a symbolic manipulation package on a computer. Taking the 
imaginary part of 5IP", the only non- vanishing contribution after performing the static 
limit comes from the singularities of the propagator. These, however, give just the delta- 
functions considered in the main text, c.f. Eq. ffTOl) . 

In the small-£ limit, the integrals may then be done analytically, and one finds 



Jim h^,k) = Z^^( 2 _^)^( 2 _l log (3 + 2V2) 

Jim Im i 7 ( u , k) = ^^^ ^yf (2 - \ log (3 + 2V2)) . (A3) 

Therefore, there is a gauge- invariant, non- vanishing 0(g 3 ) contribution to the imaginary 
part of the structure functions in the static limit. Since I find it highly improbable that 
the corresponding contribution from the second diagram in Fig. E](b) or anisotropic vertex 
corrections cancel this contribution exactly, I will refrain from evaluating also these contri- 
butions, and instead take Eqns. (IA3I) as indication that the whole 0(g 3 ) contribution does 
not vanish. 

However, both diagrams will have to be recalculated at non-vanishing frequency when 
aiming for the full LO contribution to the fluctuation coefficients /c*. 

Note added: In thermal equilibrium systems (which are necessarily isotropic) the KMS 
condition requires the imaginary part of the self-energy to vanish in the static limit (36[. Out 
of equilibrium, the imaginary part of the self-energy has to be an odd function of frequency 
[37| . but this does not seem to preclude a discontinuity at vanishing frequency. To show 
unambiguously whether there is a non-vanishing imaginary part of the self-energy in the 
static limit for non-equilibrium system, one has to calculate the above diagrams with both 
resummed propagators and vertices, which is interesting and doable, but beyond the scope 
of this work. 



[1] U. W. Heinz, '"RHIC serves the perfect fluid' - Hydrodynamic flow of the QGP," 
|arXiv:nucl-th/051205l 

[2] R. Baier, A. H. Mueller, D. Schiff and D. T. Son, "'Bottom-up' thermalization in heavy ion 

collisions," Phys. Lett. B 502 (2001) 51 [arXiv:hep-ph/000923 7 1 . 
[3] Z. Xu and C. Greiner, "Thermalization of gluons in ultrarelativistic heavy ion collisions 

by including three-body interactions in a parton cascade," Phys. Rev. C 71 (2005) 064901 

[arXiv:hep-ph/0406278] . 



13 



S. Mrowczynski, "Plasma instability at the initial stage of ultrarelativistic heavy ion colli- 
sions," Phys. Lett. B 314 (1993) 118. 

S. Mrowczynski, "Instabilities driven equilibration of the quark-gluon plasma," Acta Phys. 
Polon. B 37 (2006) 427 [arXiv:hep-ph/0511052] . 

P. Arnold, J. Lenaghan, G. D. Moore and L. G. Yaffe, "Apparent thermalization due 
to plasma instabilities in quark gluon plasma," Phys. Rev. Lett. 94 (2005) 072302 
[arXiv:nucl-th/0409068] . 

A. Rebhan, P. Romatschke and M. Strickland, "Hard-loop dynamics of non-Abelian plasma 

instabilities," Phys. Rev. Lett. 94 (2005) 102303 | arXiv:hep-ph/0412016] . 

P. Romatschke and R. Venugopalan, "Collective non-Abelian instabilities in a melting color 



glass condensate," Phys. Rev. Lett. 96 (2006) 062302 arXiv:hep-ph/0510121] . 

P. Romatschke and R. Venugopalan, "The unstable Glasma," Phys. Rev. D 74 (2006) 045011 

[arXiv:hep-ph/0605045l . 

P. Arnold and G. D. Moore, "QCD plasma instabilities: The nonabelian cascade," Phys. Rev. 
D 73 (2006) 025006 [arXiv:hep-ph/0509206] . 

A. Dumitru, Y. Nara and M. Strickland, "Ultraviolet avalanche in anisotropic non-Abelian 
plasmas," |arXiv:hep-ph/ 0604149. 

P. Arnold and G. D. Moore, "The turbulent spectrum created by non-Abelian plasma insta- 
bilities," Phys. Rev. D 73 (2006) 025013 | arXiv:hep-ph /0509226] . 

A. H. Mueller, A. I. Shoshi and S. M. H. Wong, Nucl. Phys. B 760 (2007) 145 
[arXiv:hep-ph/0607136] . 

P. Romatschke and A. Rebhan, "Plasma instabilities in an anisotropically expanding geome- 
try," Phys. Rev. Lett. 97 (2006) 252301 [arXiv:hep-ph /0605064] . 

G. D. Moore and D. Teaney, "How much do heavy quarks thermalize in a heavy ion collision?," 
Phys. Rev. C 71 (2005) 064904 |arXiv:hep-ph/0412346| . 



J. D. Bjorken, "Energy Loss Of Energetic Partons In Quark - Gluon Plasma: Possible Extinc- 
tion Of High P(T) Jets In Hadron - Hadron Collisions," FERMILAB-PUB-82-059-THY 
B. Svetitsky, "Diffusion Of Charmed Quark In The Quark - Gluon Plasma," Phys. Rev. D 37 
(1988) 2484. 

M. H. Thoma and M. Gyulassy, "Quark Damping And Energy Loss In The High Temperature 
QCD," Nucl. Phys. B 351 (1991) 491. 

E. Braaten and M. H. Thoma, "Energy loss of a heavy quark in the quark - gluon plasma," 
Phys. Rev. D 44 (1991) 2625. 

P. Romatschke and M. Strickland, "Collisional energy loss of a heavy quark in an anisotropic 
quark-gluon plasma," Phys. Rev. D 71 (2005) 125008 |arXiv:hep-ph/0408275| . 
M. G. Mustafa, "Energy loss of charm quarks in the quark-gluon plasma: Collisional vs 
radiative," Phys. Rev. C 72 (2005) 014905 [arXiv:hep-ph/0412402] . 

M. Djordjevic, "Collisional energy loss in a finite size QCD matter," Phys. Rev. C 74 (2006) 
064907 [arXiv:nucl-th/0603066] . 

S. Peigne, P. B. Gossiaux and T. Gousset, "Retardation effect for collisional energy loss of 
hard partons produced in a QGP," JHEP 0604, 011 (2006) |arXiv:hep-ph/0509185| . 
A. Peshier, "The QCD collisional energy loss revised," Phys. Rev. Lett. 97 (2006) 212301 
[arXiv:hep-ph/0605294T . 

P. Arnold, G. D. Moore and L. G. Yaffe, "Effective kinetic theory for high temperature gauge 
theories," JHEP 0301 (2003) 030 |arXiv:hep-ph/020935 3|. 

P. Romatschke and M. Strickland, "Collective modes of an anisotropic quark gluon plasma," 



14 



Phys. Rev. D 68 (2003) 036004 [arXiv:hep-ph/0304092] . 

P. Romatschke and M. Strickland, "Energy loss of a heavy fermion in an anisotropic QED 
plasma," Phys. Rev. D 69 (2004) 065005 |arXiv:hep-ph/0309093| . 

P. Jacobs, "Jets in nuclear collisions: Status and perspective," Eur. Phys. J. C 43 (2005) 467 
[arXiv:nucl-ex/0503022l . 

A. Muronga, "Second order dissipative fluid dynamics for ultra-relativistic nuclear collisions," 
Phys. Rev. Lett. 88 (2002) 062302 [Erratum-ibid. 89 (2002) 159901] | arXiv:nucl-th/0104064] . 
U. W. Heinz, H. Song and A. K. Chaudhuri, "Dissipative hydrodynamics for viscous relativistic 
fluids," Phys. Rev. C 73 (2006) 034904 [arXiv:nucl-th/0510014] . 

R. Baier, P. Romatschke and U. A. Wiedemann, "Dissipative hydrodynamics and heavy ion 



collisions," Phys. Rev. C 73 (2006) 064903 [arXiv:hep-ph/ 0602249|. 



N. Armesto, C. A. Salgado and U. A. Wiedemann, "Measuring The Collective Flow With 
Jets," Phys. Rev. Lett. 93 (2004) 242301 [arXiv:hep-ph/0405301| . 

E. Braaten and R. D. Pisarski, "Calculation of the gluon damping rate in hot QCD," Phys. 
Rev. D 42 (1990) 2156. 

A. K. Rebhan, "The NonAbelian Debye mass at next-to-leading order," Phys. Rev. D 48 
(1993) 3967 [arXiv:hep-ph/9308232| . 

J. P. Blaizot, E. Iancu and A. Rebhan, "Approximately self-consistent resummations for the 
thermodynamics of the quark-gluon plasma. I: Entropy and density," Phys. Rev. D 63 (2001) 



065003 arXiv:hep-ph/0005003|. 



M. Le Bellac, "Thermal Field Theory", Cambridge University Press 1996, sec. 2.3. 

K. c. Chou, Z. b. Su, B. 1. Hao and L. Yu, "Equilibrium And Nonequilibrium Formalisms 

Made Unified," Phys. Rept. 118 (1985) 1. 



